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Abstract. New fifth forces have not yet been detected in the laboratory or in the 
solar system, hence it is typically difficult to introduce new light scalar fields that 
would mediate such forces. In recent years it has been shown that a number of non- 
linear scalar field theories allow for a dynamical mechanism, such as the Vainshtein and 
chameleon ones, that suppresses the strength of the scalar fifth force in experimental 
environments. This is known as screening, however it is unclear how common screen- 
ing is within non-linear scalar field theories, k-essence models are commonly studied 
examples of non-linear models, with DBI as the best motivated example, and so we 
ask whether these non-linearities are able to screen a scalar fifth force. We find that 
a Vainshtein-like screening mechanism exists for such models although with limited 
applicability. For instance, we cannot find a screening mechanism for DBI models. 
On the other hand, we construct a large class of k-essence models which lead to the 
acceleration of the Universe in the recent past for which the fifth force mediated by 
the scalar can be screened. 



Contents 



1 Introduction 1 

2 Known screening mechanisms 3 

2.1 A linear theory 4 

2.2 The chameleon 5 

2.3 The symmetron 6 

2.4 The Galileon and the Vainshtein mechanism 6 

3 k-essence models 7 

3.1 k-essence screening 7 

3.2 k-essence Cosmology 8 

4 Chameleon-like screening 9 

5 Vainshtein- like Screening 11 

5.0.1 An example model with screening 12 

5.0.2 Can operators with non-trivial dependence be allowed? 14 

5.1 Multiple Operators 15 
5.1.1 The DBI case 16 

5.2 Non-perturbative treatment 17 

5.2.1 Non-perturbative examples 18 

5.2.2 Cosmology of the non-perturbative models 20 

6 Conclusions 23 



1 Introduction 

The study of non-linear scalar field theories has rapidly expanded in recent years, 
because the non-linearities enable the scalar field to be compatible with observations 
when related linear fields theories would be excluded 1 . Particular success has been had 
in constructing theories which are compatible with the null results of searches for fifth 
forces. Experimental searches for fifth forces have been conducted in the laboratory, 
and the leading results come from the Eot-Wash group [1]. The presence of fifth 
forces can also be constrained by studying the motion of bodies in the solar system [2]. 
Neither of these extremely precise searches have detected any evidence for additional 
forces, yet if a new scalar field is introduced into a theory we expect it to mediate a 
new force. 

The simplest kind of scalar field we could imagine introducing is a canonical scalar 
with a mass term. In our terminology this is a linear theory. If such a field is light and 

1 In this paper we call "linear" any model where scalar fields have equations of motion in which 
terms depend on only one power of the scalar field, all other theories will be called "non-linear" . 
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couples to matter then the fifth force experiments we mentioned above constrain the 
strength of this interaction to be at least five orders of magnitude weaker than gravity. 
This implies that the energy scale controlling the strength of the interaction must be 
five orders of magnitude above the Planck scale. This is clearly very difficult to justify 
in the context of a sensible effective field theory. 

There are a number of ways we could attempt to avoid this problem. We could 
insist that new scalar fields must be sufficiently heavy that they cannot mediate forces 
over the distance scales probed in these experiments. However attempts to explain the 
late time acceleration of the expansion of the universe typically require very light scalar 
degrees of freedom. We could insist that these new light scalar fields are forbidden from 
interacting with matter fields, for example by imposing a shift symmetry on the scalar 
field that would forbid the problematic coupling. However when studying such axionic 
models it is necessary to break the shift symmetry in order to reproduce the desired 
cosmology, and depending on how this breaking occurs the problematic coupling to 
matter may be reintroduced [3, 4]. The most successful attempts to introduce new 
light scalar fields which couple to matter without violating the results of fifth force 
experiments rely on making the theory non-linear. 

Non-linear terms can be introduced into the potential for the scalar field [5], the 
way it couples to matter [6], or the kinetic terms [7]. If the non-linearities are irrelevant 
in vacuum, but become important in the presence of sufficient quantities of matter then 
they change the way the field interacts with that matter distribution. This can result in 
the force being suppressed in experimental environments without the need to fine tune 
any parameters. Considering these theories from a classical point of view, the role of the 
non-linearities is to impede the growth of the scalar potential well around the massive 
source. If the scalar potential well is shallower than the gravitational well the scalar 
force will be weaker than the gravitational force. From a particle theory perspective 
it is possible to see why the scalar force is suppressed by perturbing the scalar field 
around a background configuration <fr . Then to second order, the Lagrangian for the 
fluctuations 5<p has the following form 



where the final term describes the interactions between the scalar field fluctuations and 
the trace of the energy momentum tensor of matter fluctuations. If the background 
scalar field configuration <po is non-trivial due to the presence of a source then the 
non-linearities present in the theory can act to make the coefficient of the kinetic term, 
Z, the mass of the field, m, or the strength of its coupling to matter, /3, vary between 
different points in space and the scalar force may be screened. If Z(4>q) increases the 
propagation of scalar fluctuations is impeded, thus making the scalar force difficult 
to transmit, this is the essence of the Vainshtein screening mechanism [8, 9] that is 
employed by Galileon models with non- linear kinetic terms [7]. If m(4>o) increases the 
scalar force only propagates over shorter distances, this is the essence of the chameleon 
screening mechanism [5]. If (3(4>o) decreases the interaction between the scalar field 
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and matter weakens, this is the essence of the symmetron screening mechanism [6] or 
the Damour-Polyakov phenomenon for dilaton models [10-12]. 

It is currently unclear whether most non-linear theories contain screening mech- 
anisms, or whether the models that have been studied are unusual exceptions. The 
non-linearities could cause the functions in the Lagrangian of Equation (1.1) to change 
in the opposite way to that required, thus increasing the strength of the force, or the 
non-linearities could cause changes in more than one function that conspire to cancel 
out the effects of screening. In this paper we study a well motivated class of non-linear 
theories and ask whether the presence of these non-linearities allows the scalar force 
to be screened. These are known as k-essence models. The k-essence Lagrangian is an 
arbitrary function of and X = (1/2) (d<f)) 2 . These Lagrangians were first introduced 
in the context of kinetically driven inflation [13, 14], and were later applied to dark 
energy [15-17]. Despite the presence of Lagrangian terms with many derivatives k- 
essence models are free of Ostrogradski's ghost and can be considered as a specific case 
of the most general ghost free scalar field theory first identified by Horndeski [18, 19]. 
A particular example is the DBI Lagrangian [20]; 



This is derived from string theory and describes the open string dynamics of D-branes. 
The DBI Lagrangian is believed to be UV complete and is therefore taken as more phys- 
ically motivated than other k-essence theories. A proof that Vainshtein-like screening 
solutions can exist for generalised Galileons, which includes k-essence models, was 
given in [21]. In this paper we specialise to k-essence and ask specifically what form a 
k-essence model must take in order to have a successful screening mechanism. 

In the following Section we will review the screening mechanisms that are currently 
known. Then in Section 3 we will introduce k-essence models and their cosmology and 
describe the approach we will take to studying screening in k-essence models. In 
Section 4 we study whether a chameleon-like screening mechanism exists for k-essence. 
In Section 5 we study whether a Vainshtein-like screening mechanism exists for k- 
essence. We give an example of a model that screens in this way in Section 5.0.1, and 
we are able to extend our analysis to the case when the k-essence Lagrangian contains 
more than one operator in Section 5.1 and then to a fully non-perturbative treatment 
in Section 5.2. In Section 5.2.2 we show that these models can also give rise to a 
cosmology that matches the one we observe. We conclude in Section 6. 

2 Known screening mechanisms 

We will take a classical approach to studying whether the scalar force is screened, 
by computing the scalar field profile around a massive source and from this deducing 
the scalar force. We will assume that quantum corrections can be safely neglected. 
In this section we review the screening mechanisms that are currently known from 
this standpoint. An important practical distinction between models will be whether 
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the non-linearities become important only inside a massive source, as in the first two 
examples we discuss, or whether they become important exterior to the source as in 
the last example. 

In this article we will only consider theories in which the coupling between the 
scalar field and matter fields takes the form 

£DA(0)T, (2.1) 

where T is the trace of the energy momentum tensor of the matter fields 2 . These scalar 
fields are referred to as having conformal couplings to matter and scalars coupled in 
this way are the most difficult to reconcile with experimental searches for fifth forces. 
The force mediated by a scalar field that couples in this way is 

J V0 , (2.2) 



M P 
where 

m = M P ^M . (2.3) 

For static, spherically symmetric scenarios the scalar force will only be non-trivial in 
a radial direction. 

2.1 A linear theory 

We begin by showing what happens in the linear case, so that it is easier to see how 
the non-linear models with screening mechanisms differ. The linear Lagrangian is 

/^-WJ'-^V + ^T. (2.4) 

The simplest scenario in which to study scalar forces is to consider a spherical, static 
massive source of density p and radius R, corresponding to a mass M c = (4/3)nR 3 p. 
Here and in what follows we assume that the background spacetime is flat and that 
the gravitational effects of the source can be well described with a Newtonian theory. 
Then inside the source the scalar field grows as 



f3p f sinh mr 



m?Mp \ mr 



(2.5) 



where, because we are interested in light scalar fields, we have made the approximation 
mR <C 1. Outside the source 

' m ' . (2.6) 



4vrM p r 



2 Other types of coupling are possible [22], and also lead to interesting behaviour, however in this 
article we restrict ourselves to the most commonly studied form of coupling. 
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Therefore outside the source the ratio of the scalar to Newtonian forces is 

= 2/3 2 (l +mr)e- mr . (2.7) 

From this is is clear that on distance scales shorter than the Compton wavelength of 
the scalar field, mr <C 1, the scalar force can only be suppressed by tuning (3 to be 
small. 




2.2 The chameleon 

The archetypal chameleon model starts from the Lagrangian 

£ , = _l( W -A'g) + gT, (2.8) 

although many other choices of potential are possible. The role of the potential is to 
make the effective mass of the scalar field much larger inside a massive non-relativistic 
source of density p than outside, this is because the scalar field dynamics are governed 
by an effective potential 

The mass of the field at the minimum of the potential is 

"'(a^F' < 210 > 

Clearly the mass of the field changes with the local value of p. 

How to solve the chameleon equations of motion has been described in detail 
in [5, 23], but the important features of the solution are that inside a sufficiently 
massive object the scalar field takes a constant value everywhere except very close 
to the surface. This is how we see that the scalar potential well is shallower for the 
chameleon, than for a corresponding linear scalar field, because the growing mass of 
the scalar field makes it increasingly costly for the scalar field to decrease its value. 
Outside the massive source, and at distances smaller than the Compton wavelength of 
the field rm^ <C 1, the ratio of the chameleon scalar force to the Newtonian one is 

n _ ,9(000 - C ) 8nM 2 P R 

F N ~ M P M c ' ^ 

where C = ^A 5 /M P /3p and 0^ = y/A 5 /M P /3p 00 are the values of the scalar field that 
minimise the effective potential inside and outside the object respectively. From this 
expression it is clear that the scalar force is weaker than gravity if the depth of the 
scalar potential well inside the source, (3A(fi/M P = f3(cj) QC> — <fi c )/M P , is shallower that 
the corresponding gravitational potential well, $ = M c /8nMpR. 
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2.3 The symmetron 

The original symmetron model begins with the Lagrangian 

C, = -\{d^ + ^-X^f-^- 2 T. (2.12) 

As with the chameleon model the position of the minimum of the effective scalar 
potential varies depending on the local energy density. When this energy density is 
sufficiently high p > M 2 p 2 the minimum is at <fi = 0. At this minimum the coupling 
between scalar fluctuations and matter vanishes, and so the scalar field no-longer re- 
sponds to the presence of a massive source. In the interior of such dense objects the 
scalar potential well stops growing, and so the scalar force sourced by that object is 
suppressed. Outside the same source object of density p and radius R that we consid- 
ered before, and inside the Compton wavelength of the field pr <C 1, the ratio of the 
symmetron scalar force to the Newtonian one is 

p 2 8ttM 2 R ( R\ 

At infinity the minimum of the effective potential for the scalar field is at = 
p/\^X and within a sufficiently massive object the minimum is at <p c = 0, so as in the 
chameleon case we see that the scalar force is suppressed compared to gravity whenever 
the potential and form of the coupling to matter conspire to stop the scalar potential 
well, A(j) 2 /M 2 , from growing within the source object as fast as the gravitational well 
grows, $ = M c /%ixM 2 P R. 

2.4 The Galileon and the Vainshtein mechanism 

The final screening mechanism that is known behaves somewhat differently from the 
two previously discussed. The simplest Lagrangian that displays this form of screening 
is 

^=->) ! -^(W + ^. (2.14) 

In this example the non-linearities are gradient terms, and therefore they can become 
important in the exterior of the source object. The Galileon theory shows how such 
higher order derivative interactions can be constructed without giving rise to ghosts. 
The non-linearities dominate the behaviour of the scalar field within a specific radius, 
known as the Vainshtein radius; 



(2 , 5 ) 



1/3 

\2ttM p J A 

On solving the equations of motion we find that within this radius the ratio of the 
scalar to Newtonian forces is 



2^ — ] (2.16) 



- 6 - 



Outside the Vainshtein radius the scalar force is unscreened, but within it the force is 
suppressed, and the further inside the radius the more the force is suppressed. 

In the discussion of screening mechanisms for k-essence models that follows we 
will consider two possibilities, one that the non-linearities may only be important 
inside the source, we will call these chameleon-like screening mechanisms, and second 
that the non-linearities can be important outside the source object, we will call these 
Vainshtein-like screening mechanisms. 



3 k-essence models 



We will consider an operator expansion of the k-essence Lagrangian, C(<f), X). It may be 
that to describe a particular theory completely requires an infinite series of operators, 
however in any given situation we expect only a finite set of operators to dominate 
the behaviour of the field, and we are interested in which of these operators allow 
for screening. We assume that the scalar field couples to matter conformally. Such 
a coupling is always present unless there is a symmetry to forbid it, and this is not 
the case for a general k-essence scenario. Therefore our starting point is the scalar 
Lagrangian 

= £ <7>o 4m ^ m - ^ T , (3-1) 

n,m ^ 

where 

X = ^) 2 , (3.2) 

and T is the trace of the energy momentum tensor for matter. We assume that the 
coupling to matter and the non-linearities of the theory are controlled by different mass 
scales; Mp/j3 and M n ^ m respectively. 

The simplest scenario in which to study the existence of absence of screening is 
a spherically symmetric, static non-relativistic source of density p, radius R and mass 
M c . Therefore T = —p(r) = —pQ(R — r). Working around a flat space background 
the resulting equation of motion is 



m—1 



(3.3) 

where a/— g is not automatically trivial, for instance when spherical coordinates are 
used. 

3.1 k-essence screening 

Our approach in what follows will be based on the examples of known screening mech- 
anisms discussed in the previous section. We will assume that there is a near vac- 
uum region in which the non-linearities of the theory are unimportant. If a massive 
source object is introduced into this vacuum then we expect the non-linear terms of 
the theory to be excited in the neighborhood of the source, and provide a dynamical 
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screening mechanism. As mentioned before we will consider two cases separately one 
where the non-linear terms are excited outside the source object, which we will call 
Vainshtein-like, and then other where they are only excited inside, which we will call 
chameleon-like. 

We will solve the equations of motion in a piecewise fashion, breaking space 
into regions inside and outside the object and regions where the non-linearities are 
important and un-important. We will impose that the scalar field and its derivative 
are continuous when passing between these regions. In what follows we also assume 
that there is only one region of space where the non-linear terms are important, and 
that within this region there is a power-law solution to the equation of motion. It is 
clear that these last two assumptions are a considerable simplification of what may be 
possible, we are clearly only studying the presence or absence of the simplest types of 
screening mechanism in k-essence theories. We leave a generalisation of this to more 
complicated cases for future work. 

There is a very clear distinction between theories where the non-linearities are 
purely in the potential or in the form of the coupling to matter, and those where the 
non-linearities rely on derivative operators. When the non-linearities depend only on 
</> and not on its derivatives then the equations of motion will have the form 

□0 = /(<M, (3-4) 

where the function / depends on the theory being considered. Examples of such the- 
ories include the chameleon and symmetron models discussed in the previous section, 
and in those cases f(<fr,p) can always be locally approximated by a (potentially tachy- 
onic) mass term. This makes the equation of motion homogeneous, and the solutions 
have particularly nice properties and contain two constants of integration which are 
fixed by imposing boundary conditions. In this article we focus on terms with a non- 
linear dependence on X, and their role in the screening of a scalar force, and therefore 
we will not be able to make a homogeneous approximation to the equations of motion. 
The power law approximation that we choose to make always allows us to find one 
solution to the equations of motion, but there may be others that we miss with this 
approach. Systems described by equations of motions of the type given in Equation 
(3.4) have been well studied in the literature, hence in this article we focus on sys- 
tems with kinetic non-linearities in X. This means that the following analysis will not 
include the chameleon and symmetron cases. 

Clearly many of the operators in Equation (3.1) are non-renormalisable. In this 
paper we take a bottom up approach to identifying screening mechanisms, asking 
only; which theories are able to screen the classical scalar force? If we do find that 
such screening mechanisms are possible then it will require further study to identify 
whether they can be embedded in a well behaved effective field theory. 

3.2 k-essence Cosmology 

We digress briefly to discuss the role of k-essence fields in cosmology, k-essence models 
have been studied as explanations for the late-time acceleration of the expansion of 
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the universe. Commonly studied models have a separable form 



C D V^K(<p)p(X) . (3.5) 
If the scalar field is treated as a perfect fluid, then it has equation of state 

= 7 r V J^ , (3.6) 

where p } x indicates a derivative of p with respect to X. Clearly with a suitably chosen 
p the scalar can have an appropriate equation of state to allow it to act as dark energy. 
When Ucf, is sufficiently close to —1 the energy density in the scalar field will redshift 
away more slowly than the energy density in matter. The instant when the scalar field 
comes to dominate the dynamics of the universe is controlled by the energy scales of 
the higher order operators in the Lagrangian, and to solve the coincidence problem 
it is necessary to tune these energy scales to take values corresponding to the energy 
density of the universe today. We will give a specific class of examples of this in Section 
5.2.2. 



4 Chameleon- like screening 

Chameleon-like screening happens when the non-linearities of the theory are only im- 
portant inside the source. To study this we consider scenarios in which only one type 
of non-linearity dominates. If we assume that the density of the source is a constant 
p inside a radius R and that outside there is vacuum, then for r < R the equation of 
motion becomes 

r BM n+4m ~ 4 or 3 
n / dr r 2 ct) n - l X m + - P — = mr 2 0V* . (4.1) 



3M F 

We choose a power law solution 

= Cr a . (4.2) 

Then matching the powers of r in the equation of motion requires 

2m . , , 

a = , 4.3 

2m + n - 1 v ' 

and a full solution of the equation of motion requires 

C 2m + n-l = 3 X 2™- 2 pM n+4m ~ 4 M c 

ira 2m (6m + 2n-3)M P R 3 ' 1 ' ' 

For the force 0' to have a smooth behaviour at the origin we need to impose a > 1. 
This requires 

— |m| < n + m — 1 < \m\ . (4-5) 
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As we move further away from the origin the non-linear terms become less im- 
portant and we recover the linear behaviour. We will call the radius at which we cross 
from the non-linear regime to the linear one R$ and it must satisfy 



M 



\r=R S 



1 . 



(4.6) 



Evaluated on the solution of Equation (4.2), Equation (4.6) determines Rs in terms of 
the Lagrangian parameters 



CR a s ~ 2 



(3p 



a 2 (Qm + 2n-3) M p 



(4.7) 



Whilst still in the interior of the source object, but outside the region where 
non-linear terms dominate, the solution is 



Ppr 2 , a 
bMp r 



and the exterior solution is 



+ 



d 



(4.8) 



(4.9) 



where a, b, d, (p^ are constants, to be determined by boundary conditions and continuity 
of <f) and its derivative at R and Rs- The constant d controls the strength of the force 
exterior to the object. This is weaker than gravity when 



(3d 



Solving the continuity equations we find 

(3p 



< 



Mr. 



8ttMJ, 



(4.10) 



d 



(Rs - R s 



CaR a s +l 



3M P 

Substituting in the form of C from Equation (4.4) and defining a new length scale 

6 



(4.11) 



E>3 p3 



a(Qm + In — 3) 



(4.12) 



this can be more concisely written as 

/3M C 



d 



5- 



AttMp \R 3 

and therefore the scalar force is weaker than gravity when 



(4.13) 



R s 
i? 3 



1 



< 1 . 



(4.14) 
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This is the k-essence version of the chameleon thin-shell condition. It is clear that 
when n, m ~ 0(1) the thin shell condition in Equation (4.14) can only be satisfied by 
objects with a specific correlation between their mass and radius. This cannot include 
all objects used in tests of gravity, which cover a wide range of densities and radii. 
Therefore the thin shell condition in Equation (4.14) tells us that we cannot screen 
the fifth force around a wide variety of objects unless m and n are very large. This 
requires that inside the source object an operator that we would normally consider to 
be highly irrelevant dominates the dynamics of the scalar field. It is difficult to see 
how such a theory could be constructed in a well controlled way. 



5 Vainshtein-like Screening 

We now consider the Vainshtein-like situation where the non-linearities of the theory 
are important exterior to the source object. We can treat the massive source as a delta 
function p(r) = M c 5(r). Assuming that only one of the non-linear terms dominates 
close to the matter source the equation of motion becomes 



/3M C 
4vrM ; 



JO 



Clearly the leading contribution from the right hand side of this equation must be 
independent of r in order to balance the source on the left hand side. Therefore, away 
from the origin at r = 0, the two terms on the right hand side of Equation (5.1) 
must have equal derivatives with opposite sign. As we are interested in power law 
solutions for 0, the integral on the right hand side of Equation (5.1), cannot give a 
constant contribution that is independent of r. Therefore the leading order solution to 
Equation (5.1) has to have n = 0. If we were considering a DBI theory the condition 
n = would mean that the tension function T(0) would have to remain constant over 
the range of <fi necessary to describe the configuration around a massive source. 

After setting n = we need the right-hand side of Equation (5.1) to be a constant. 
Assuming that is well approximated by a power law 



then this occurs when 



<j> = Cr a , (5.2) 

2m — 3 , 

a = . 5.3 

2m - 1 v ' 

The equation of motion (5.1) is then satisfied if 

aC = ^— . 5.4 

For this to be a consistent solution to the equation of motion we need the non- 
linear operator we are interested in to dominate over the canonical kinetic term in the 
equations motion. This requires the following consistency condition 

1< mM^-^lX™- 1 ! . (5.5) 
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For the power law solution this requires 

m— 1 



rim - m - a) < R i r m- a) = ml/2 m " (5 6) 



The radial force mediated by this scalar field is 

* ~ iWp ~~ Wp 



(5.7) 



For the moment we do not impose that the force does not diverge as r — > because we 
have treated the source as a point-like object and the behaviour of the field as r — > 
will be modified when the finite extent of the source is taken into account. We will 
return to this in Section 5.2. This k-essence force is weaker than gravity when 

r 8nM P p\aC\ ' ^ 

For a successful screening mechanism the force should be weaker than gravity close to 
the source at small r. Therefore we need a + 1 > 0. This requires 

m > 1 or m < — . (5-9) 

If we are only interested in integer values of m, then this makes no restriction. 
So we have two length scales that control the behaviour of the field 

Rg= \¥^) {snMpM* ) ' (5 - 10) 

within which the scalar force is weaker than gravity, and 

where the non-linear terms dominate over the canonical ones. A successful screening 
mechanism must have Rg < Re, which requires 

(2/3 2 ) 2 ™- 1 > 1 . (5.12) 

5.0.1 An example model with screening 

To illustrate more clearly how screening occurs for k-essence theories we choose that 
the non-linear term that becomes dominant has m — 2. The theory we consider is 

C=-X*%-„. (5.13) 

M 4 M P H v ; 

The resulting equation of motion is 

n 0+^W) 2 M-^ = O. (5.14) 
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For a static, spherically symmetric configuration sourced by a point like object of mass 
M c at the origin, we can integrate the equation of motion once to find 



4>' + -A' 3 = ^- 2 . (5.15) 

If the non-linear term dominates close to the source then this has solution 

= CV 1/3 , (5.16) 

where 

C = 330M ' M ' ■ (5.17) 
The two distance scales we are interested in are 

^ 2l/ ^ /2 (wf- (5 ' 19) 

and the consistency condition that Rq > Rq requires 

2(3 2 > 1 . (5.20) 
The ratio of the k-essence scalar force to the Newtonian force is 

7-i / \ 4/3 / \ 4/3 

Fa, r \ 1 ^^,if T 



2^ - (5.21) 

Therefore the scalar force becomes more screened as we travel further within the radius 
Rq. This behaviour is very similar to that of the Vainshtein case considered in Section 
2.4, but k-essence screening occurs more slowly with r than the Vainshtein example 
discussed in Section 2.4 and the screening radius has a different dependence on the 
couplings of the scalar field and the mass of the source. 

One system within which the scalar force must be screened is the interactions 
between the Earth and the Moon. We have very good experimental evidence that the 
orbit of the Moon around the Earth is well described by General Relativity. Therefore 
if M c = M ffi then Rq must be bigger than the distance from the Earth to the Moon, 
approximately 3 x 10 5 km. This imposes 

PM < eV . (5.22) 

This is a surprisingly low energy scale. It is intriguing that this allows for M ~ 10 -3 eV, 
the dark energy scale, suggesting that the higher order operators could be related to the 
currently unknown mechanism which controls the value of the cosmological constant. 
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5.0.2 Can operators with non-trivial <j) dependence be allowed? 

In order to satisfy equation (5.1) we had to choose n = 0. This seems rather restrictive, 
and we now ask whether the scalar field is allowed to enter into the Lagrangian without 
derivatives if two (or more) non-linear terms become important at the same time and 
conspire to cancel the unwanted integral terms in the equation of motion (5.1). The 
equation of motion is now 



4vrM 



p 



- j dr r 2 (M 4 - n - im n<p n ~ l X m + N 4 - q ~ 4p q ( f) q ~ 1 X p ) (5.23) 
+ [r 2 d r <f){mM 4 - n - 4m <f) n X m - 1 +pN i - q - ip <p q X p ~ 1 )}* , 



where we have allowed the energy scales controlling the two nonlinear terms, M and 
N to be distinct. 

We still need the integral on the right hand side of Equation (5.24) to vanish, 
which means the solution must satisfy 

6 n -i = — X p ~ m . (5.24) 

Then satisfying the equation of motion for a power law solution <fi = Cr a requires us 
to choose 

= 2 P~ 3 (5.25) 
q + 2p - 1 ( J 

and 

cq+ 2 P -i = 2 p -*n(3M c 

7ia 2p - l (pn - mq)M p N 4 -i- 4 P V ; 

The condition for the vanishing of the integral part of the equation of motion, 
given in Equation (5.24) is satisfied on these solutions if 

Qn-q-2(p-m) r a{n-q)-2(a-l){p-m) H 2 m_ P Q; 2 (P _m ) 27) 

nM 4-n-4m ' 

and the solution to this must be independent of r, therefore 

a(n-q) = 2(a - l)(p-m) . (5.28) 

We also want the relationship between M and iV to be independent of the mass of 
the source, in order to ensure that the equation of motion can be satisfied for a vari- 
ety of different source objects. This means that the consistency condition has to be 
independent of C. This requires 

n-q = 2(p-m). (5.29) 

Unfortunately satisfying equations (5.28) and (5.29) is not possible without choosing 
p = m, which would mean that the two operators are identical. 
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5.1 Multiple Operators 

Let us now consider the case when the number of relevant operators is finite or even 
infinite, with each only depending on X and not on 0. We define 

M 4 " 4m = c m M 4 " 4m , (5.30) 

where M is a characteristic scale and the Lagrangian operators under consideration 
are 

X m 

°m = C m M4m _ 4 ■ (5.31) 

We will analyse this situation by singling out one term defined by its index m and 
compare it to all the other operators. If this operator dominates then the solution to 
the Klein-Gordon equation is 

<j>(r) = C m r a ™ , (5.32) 

where a m and C m have already been defined in Equations (5.3) and (5.4). We define 
the radii Rp tTn to describe when the O m operator dominates over other terms in the 
equations of motion: 

O m > O p for m > p if r < R m>p , (5.33) 

and 

O m > O p for p > m if r > R ni:P . (5.34) 

The first condition involves a finite set of operators with fewer powers of X while the 
latter may involve an infinite number of conditions. An expression for R p ^ m can be 
found using the power law solution computed in Section 5 when the O m operator is 
dominant 

(2m-l)/2 

(5.35) 

For simplicity we have defined the dimensionless combination 



a 2(rn- p) = (5 36) 

mc m 

Using the expressions (5.3) and (5.4) for a m and C m respectively Equation (5.35) can 
be rewritten as 

D _ ol/4 (l-2m)/2 ( P M c \ /, o 7 X 

There is also a radius R m ,G given by Equation (5.10) within which the scalar 
force is suppressed compared to the gravitational force, this is closely related to the 
Vainshtein radius discussed in Section 2.4. The scalar force is successfully screened 
when 

< Rm,G < R m , (5.38) 

where we have defined 

i?+ = inf Rp, m , (5.39) 

p>m 
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and 

R~ = sup R Pim . (5.40) 

When this is not the case then the operator O m does not lead to a Vainshtein mech- 
anism. Notice too that for any given operator the Vainshtein mechanism does not 
extend to arbitrarily small objects: when r < i?~ the approximation breaks down 
and O m is not the leading operator anymore. Moreover, nothing prevents multiple 
operators O m from satisfying this inequality for objects with a size R < R G . In this 
case, the theory which is highly non-linear has many branches. The inequality (5.38) 
is satisfied provided 

P~ < < /3+ , (5.41) 

where 

/3- = 2- 1 / 2 supa^- 1 m 1 / 2 ( 2m - 1 ), (5.42) 

and 

/3+ = 2" 1 / 2 inf a2- l m l ^ 2m -V . (5.43) 

p>m 1 

Hence a given operator O m leads to a Vainshtein mechanism only for a limited range 
of couplings. If this last inequality is not satisfied, no Vainshtein mechanism is at play. 

5.1.1 The DBI case 

The operator expansion for the DBI model corresponds to 

As we have seen there is no Vainshtein-like screening mechanism when there is a non- 
trivial dependence on 0. Therefore the Vainshtein mechanism is only compatible with 
models where both T(4>) and V((p) are constant over the range of field values probed 
by the experiment. In this case setting M = T 1//4 we have a theory of the form given 
in Equation (5.30) 

c m = l 1 1 -( 1 --l)...( 3 -- m ) = (-1)- , r( :7> ■ (5.45) 
m!2V2 / V 2 / T(-l/2) T(m + 1) V ; 

It can be easily seen that a pm is a decreasing function of p for any m. This implies 
that 

R m = Rm-l,m , (5.46) 

and 

(5.47) 

Now we find that 

P~ m = 2-V*o£ 1 ro 1 / 2 < am - 1 > , (5.48) 
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and 

fa = 2^ 2 lim a^m 1 ^ 2 ™- 1 ) . (5.49) 



pm 

p— >oo 



Therefore 

fa > fa , (5-50) 

which contradicts the requirement of Equation (5.41) and so there is no Vainshtein-like 
screening mechanism at play for DBI theories. 

5.2 Non-perturbative treatment 

With the understanding we have gained from studying k-essence theories operator by 
operator we can now extend this to a non-perturbative study. We will also be able to 
allow the source object to have finite extent. Requiring that the solution be smooth 
inside as well as outside the source will be even more restrictive on the allowed set of 
solutions. The previous considerations have taught us that the Vainshtein mechanism 
can only exist when the leading non-linear terms in the Lagrangian have no non-trivial 
4> dependence. So we start from a Lagrangian 

C = M',( 2 4)-i±T, (5.51) 
J \M 4 J M P ' v ; 

where / is an arbitrary function. In particular, this requires that any potential terms for 
quintessence must be subdominant over the range of field values probed by gravitational 
experiments. For these theories, the equations of motion in a spherical setting read 
simply 

d ( 2 d(p A (3pQ(r — R) 2 



where /' = ^j^- This can be integrated to give 




dr 



' 2 \ BM r I -m, r < R , 

(5.53) 



8ttM p ii r>R 



where we have assumed that (d<ft/dr)f' is continuous at r = R and that the solution 
does not diverge as r — )• 0. 

It is easier to study solutions of this equation in terms of dimensionless variables, 
therefore we define 

x = L , (5.54) 
<p = M 2 Ru(x) . (5.55) 



and 

We also define the distance scale 
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which plays an important role in determining the dynamics of the scalar field solution. 
For any continuous /, one can invert (5.53) and find a function g(y) that satisfies 



g(y)f'(g (v)) = v 



(5.57) 



where g = du/dx and y = (R±/R) 2 x. Notice that g is an odd function of y and can 
always be written 

g(y) = y~g(y 2 ) • (5.58) 

This will be useful later when we discuss the cosmology of these theories. Once this 
inverse has been identified we can solve for the profile inside the object, r < R 



u = u(0) 



R+ 
R 



-i 2 



x dx 



(5.59) 



Outside the object, r > R, the profile is 



u = u(0)+ / g 



R* 
R 



x dx 



R* 
~R 



-i 2 



.f 2 



dx 



(5.60) 



Outside the object the scalar force is weaker than the gravitational force provided 



du < j_(R 1 
dx ~ (3 2 \R J x 2 



or equivalently 



9 



R* 
R 



x 2 



1 (R 



P 2 \ R J x 



(5.61) 



(5.62) 



5.2.1 Non-perturbative examples 

First of all, we consider the power law case discussed earlier in Section 5. If 



x 



fix) = x H 



when the non-linear terms dominate, we find 



g(y) 



2 m y 



in 



l/(2m-l) 



which gives rise to 



(r) = MR 



R 2 xl/(2m-l) (s^s)/^!) 

) \r) 



(5.63) 



(5.64) 



(5.65) 



This nicely recovers the solution we obtained previously in Equations (5.3) and (5.4). 
Equation (5.62) applied to this choice of / reproduces the expression for Rq, the radius 
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within which the scalar force is suppressed with respect to gravity, from Equation 
(5.10). From Equation (5.53) it is clear that d<fi/dr — > as r — > 0. Therefore in the 
very center of the massive object the linear terms must dominate the behaviour of 
the scalar field. However as d(p/dr grows rapidly the non-linear terms soon come to 
dominate the behaviour of the field, and continue to do so until r ~ R*. 
As a second example we take DBI with 



implying that 



f(x) = VTTx , (5.66) 



g(y) = -7=== 2 ■ (5-67) 



The solution inside the source is 



0W = M 2 ^^l- V /l-4^| . (5.68) 
When R± > R the scalar force, which is proportional to d<p/dr, will diverge at 

(R \ 2 
~\/^R~ ) ' (5.69) 

which lies inside the source object. Clearly this is not a physically acceptable situation. 
Therefore we again conclude that there is no screening mechanism for DBI models. 
If we extend this analysis to models of the form 



f(x) = Vl + x m , (5.70) 

we have now 

rug 2 ™" 1 

V ■ (5.71) 



2 y /l + g* m 

If there were singularities present as in the DBI example these would correspond to a 
divergence of g. In the case when g > 1 we have 

777 

jg™- 1 = y ■ (5.72) 

When m ^ 1 it is not possible for g to diverge at finite y, which we can equivalently 
state as there are no divergences in the scalar force at finite r. It is clear from Equation 
(5.72) why this does not apply to the DBI example with m = 1. We conclude that for 
models with m > 1 it is possible to have a successful screening mechanism. 

Finally we analyse a class of models related to DBI that will prove to have an 
interesting cosmology 
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where h(x) = J2T= 



CnX" 



From Equation (5.53) we find that d<p/dr — > as r — > 0, 



implying that deep inside a dense body 

f(x) ~ \/l + CtX. 



(5.74) 



As we have shown for the DBI case, this approximation breaks down inside the body 
as d<j)/dr increases unboundedly. When this happens the higher order terms in the 
polynomial h(x) start dominating and we can approximate 



/(*) 



m/2 



(5.75) 



This happens as soon as cp' ^> M 2 and is guaranteed (see (5.65)) to be the case as long 
as R+ 3> R and r < 0(R+). Close to R+, operators of lower order begin to dominate 
again. In this case, the screening mechanism operates for distances larger than R and 
smaller than Rq which is of the same order of magnitude as R±. In the following 
section we will show that if we want this scalar field to be connected to the late-time 
acceleration of the expansion of the Universe then M ~ 10 -3 eV. This implies that 
solar system objects are all screened in these k-essence models. 

5.2.2 Cosmology of the non-perturbative models 

We now digress to discuss the cosmology of the k-essence models that allow for a 
Vainshtein-like screening mechanism. Taking a flat FRW metric 



ds 2 



-dt 2 + a 2 dx^ 



(5.76) 



the equation of motion for a homogeneous scalar field described by Equation (5.51) 
can be expressed after integration as 



P Pr 



M 4 



2M f 



(5.77) 



where we have assumed that <ft does not diverge as the scale factor goes to zero in 
the early Universe and, because only non-relativistic matter contributes to the scalar 
equation of motion, we can assume that p m ~ a -3 . The solution to this equation can 
be obtained as the analytic continuation <j)/M = —ig(—ij3p m t/2MMp) or equivalently 



2M P 9 



2MMf 



(5.78) 



where g was defined in Equation (5.58). The Friedmann equation for H 
determines the whole cosmological evolution 



H 2 



3M 2 



P<t> 



Mr 



Pm + Pi 



then 



(5.79) 
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where p r is the energy density of all the relativistic species and the energy density of 
the scalar field is 

^f(-w)+ Ml i{-w) ■ <«■») 

We are particularly interested in the late time evolution of the Universe. This can be 
conveniently analysed using the equation of state 

w+ = ^ , (5.81) 



where 



Using the equation of motion (5.77) we find that 



M 4 /(- 



ws = ■ Jy M4J ., . (5.83) 



The scalar field dynamics mimic a cosmo logical constant ~ — 1 when £ m <C 

Let us give a class of models leading to late time acceleration. We choose to start 
with a theory of the form 

where h(x) = YLT=i c « x S sucn theories were discussed in the previous Section. We 
assume that in the recent universe, X <C M 4 , then the long time solution to the 
equation of motion is 

Y 2 Cl M P v ; 

In this case the equation of state becomes 

w * = ~ WML i 1 ' (5 - 86) 

2ciM 4 M'£ ' 

We are interested in when the scalar equation of state approaches that of a cosmological 
constant with w = —1. We define the time t$ to be when the equation of state of the 
scalar field is w = —0.99 in the recent past of the Universe, and find that it is given 
by the expression 

T ~ 7a ^mhF p ' (5 - 87) 



or equivalently 



f~ 9.5-1*. (5,8) 
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When this happens the energy density becomes p$ ~ M 4 which we want to identify 
with 3Q\HqMp in order that the scalar field explains the late time acceleration of the 
expansion of the universe. Therefore 

*+ 10/3 ^ m0 (5.89) 



where we have set to-^o ~ 1 an d ignored other order one factors. Imposing that this 
happens for a redshift of order one requires 

ci ~ 2 6 /3 2 . (5.90) 

We also need to check that the presence of this scalar field does not disrupt the 
cosmology of the early universe. For the same reasons that the higher order terms in 
the Lagrangian dominate inside the dense sources used in fifth force experiments, the 
higher order terms will also dominate in the early universe. This gives the following 
equation of motion for the scalar field at early times 

&Pmt (5.91) 



M 2 J mM 2 M P ' 

where m = 2p, c m > and p needs to be odd to guarantee that the kinetic energy in 
the Lagrangian should be positive. In what follows we will specialise to the case m = 2 
for algebraic simplicity, but the full equations are straightforward to derive. 

In the early universe the ratio of the scalar density to the matter density is 

* = P2pmt \ . (5.92) 



p m 4^M 

During the matter dominated era p m ~ t~ 2 and so the ratio of scalar and matter 
densities remains constant and less than one if 

fi 0m /3 2 < • (5-93) 

Close to the current time the higher order operators in the scalar field Lagrangian 
become smaller and no longer dominate the scalar field dynamics, then the field starts 
to behave as a dark energy candidate as discussed above. 

At even earlier times during radiation domination, the matter energy density 
scales as p m ~ t~ 3 ^ 2 therefore p m t 2 ~ t 1 / 2 is an increasing function of time. If the 
scalar field energy density is subdominant during the matter epoch then it is also 
forced to be subdominant in the radiation dominated epoch. 

On top of having a time dependent equation of state which departs from -1 in the 
past of the Universe, the effects of the scalar field on the growth of structure could also 
be interesting. Indeed for clusters of masses between M c = 1O 12 M and M c = 1O 15 M , 
the scale 

^(J^T, (5.«) 



8irM P M 2 



-22- 



which gives the order of magnitude for scales below which the fifth force is screened, 
varies from 0.5 Mpc to 15 Mpc (for (3/8ir ~ 1) implying that the effect of the scalar 
field could be felt on the largest scales. The study of this phenomenon is left for future 
work. 

6 Conclusions 

In this article we have asked whether k-essence theories possess a screening mechanism. 
We have been unable to find a chameleon-like mechanism where the non-linear terms 
are only important in the interior of a massive source, however we find that it is possible 
for k-essence models to have a Vainshtein-like screening mechanism for which the non- 
linear terms are important exterior to the source. However screening is not possible 
for all k-essence models, and in particular such a screening mechanism does not exist 
for the well motivated DBI model, k-essence models that are suitable for screening 
are dominated by operators which have no <j) dependence. For models that do screen 
we find that the non-linear screening dynamics become important on distance scales 
below 



and that within this radius the fifth force mediated by the scalar field rapidly become 
smaller than the gravitational force. 

In addition we find that there are models where the fifth forces can be screened, 
which also give rise to a cosmology that matches the one that we observe in our universe. 
The same non-linearities that screen the scalar fifth force in the presence of massive 
sources, also cause the effects of the scalar field to be suppressed in the early universe. 
At late times the equation of state can approach w = —1 and so the scalar field behaves 
like a cosmological constant term. Recovering the observed cosmology requires tuning 
one of the parameters in the scalar theory to the scale of the observed cosmological 
constant and we have nothing to add to the discussion of how such a hierarchy of 
scales could be achieved. We have shown that a k-essence model controlled by such 
low energy scales can have a successful screening mechanism to hide from fifth force 
searches and a successful background cosmology. At the perturbative level, we expect 
effects to appear in the Mpc range with possible observational consequences. This is 
left for future work. 
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